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EXECUTIVE  SUMMARY 

A  new  computational  procedure  for  including  welding-induced  residual 
stresses  and  distortions  from  weld  simulations  in  the  SYSWELD  software  code  in 
structural  Finite  Element  Analysis  (FEA)  simulations  performed  in  the  Abaqus 
FEA  code  is  presented.  The  translation  of  these  results  is  accomplished  using  a 
newly  developed  Python  script.  Full  details  of  extracting  the  results  from 
SYSWELD,  translating  these  results  using  the  newly  developed  Python  script,  and 
including  the  results  in  Abaqus  are  provided  through  an  illustrative  example 
analysis  of  a  butt-welded  aluminum  plate.  The  development  of  this  capability  is 
motivated  by  the  Navy’s  increased  use  of  lightweight  aluminum  structures  where 
including  the  effects  of  welding  in  structural  components  can  be  critical  to 
understanding  structural  response.  The  integration  of  SYSWELD  and  Abaqus 
enables  Code  60  engineers  to  increase  the  efficiency  of  design  through  the 
integrated  computational  materials  engineering  (ICME)  paradigm,  which 
emphasizes  the  use  of  computational  tools  to  understand  design  performance, 
while  minimizing  the  amount  of  traditional  physical  testing.  This  project  is  part  of 
a  program  entitled  “Integrated  Computational  Design  and  Analysis  -  Aluminum 
Ship  Materials  and  Structural  Performance  (Al  ICDA).  ” 


BACKGROUND 

As  a  result  of  increasing  requirements  for  high-speed  vessels,  the  U.S.  Navy  is  increasing 
use  of  lightweight  structures,  most  notably  in  aluminum-based  construction,  where  welding- 
induced  stresses  and  distortions  can  have  a  significant  impact  on  structural  behavior.  The  effect 
of  these  welding-induced  quantities  often  manifest  themselves  on  a  structural  level,  such  as  crack 
development  with  no  appreciable  global  structural  load  as  a  result  of  residual  stresses  or 
premature  buckling  under  compressive  loading.  In  combination  with  this,  many  of  these  high¬ 
speed  structures  exhibit  materials,  welding  configurations,  and  structural  details  that  differ  from 
historical  U.S.  Navy  structures.  As  such,  it  is  difficult  to  perform  evaluations  or  improve  the 
design  of  these  structures  using  well-developed  “rule  of  thumb”  criteria. 

Since  a  full  testing  program  examining  the  performance  of  each  detail  under  typical  naval 
shiploads  is  not  economically  feasible,  extensive  computations  accompanied  by  targeted  testing 
are  required  for  evaluation  and  development  of  new  designs.  However,  in  order  to  perform  such 
computations,  an  efficient  method  to  include  welding-induced  stresses  and  initial  distortions  in 
these  computations  is  required.  This  approach,  falling  within  the  Integrated  Computational 
Materials  Engineering  (ICME)  paradigm,  hinges  on  both  validation  and  verification  (V&V)  of 
the  applied  methods  and  the  seamless  communication  of  the  computational  tools  of  interest. 

Historically,  computational  tools  for  analyzing  welding  and  structural  behavior  have  been 
used  in  an  isolated  setting.  In  the  present  effort,  a  process  to  transfer  results  from  weld  process 
simulations,  performed  using  the  SYS  WEED  welding  simulation  Einite  Element  Analysis  (EEA) 
software  package,  to  coupon  and  structural  scale  calculations,  performed  in  the  widely-utilized 
Abaqus  EE  software  package,  is  described.  While  Abaqus  has  been  used  extensively  for 
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simulating  the  mechanical  behavior  of  aluminum  structural  components  [1-4],  no  documented 
integration  of  SYSWELD  and  Abaqus  was  identified. 

Overview  of  SYSWELD  Software 

SYSWELD  is  a  commercially- available  EEA  software  package  developed  by  ESI  Group 
for  use  in  thermo-mechanical  and  thermo-metallurgical  welding  process  simulations.  Specific 
welding  processes  can  be  directly  simulated  with  the  ability  to  input  welding  parameters, 
including  travel  speed,  arc  voltage/current,  weld  pool  size,  and  clamping  conditions.  SYSWELD 
uses  material  databases  that  contain  the  thermal,  metallurgical,  and  mechanical  material 
properties  for  the  materials  being  simulated.  Though  widely  used  and  accepted  for  capturing  the 
welding  process  in  steel,  there  are  relatively  few  examples  of  SYSWELD  used  to  simulate 
aluminum  welding  processes  [5-7].  The  reader  is  referred  to  the  SYSWELD  Toolbox  for  full 
details  on  the  usage  and  input  format  for  SYSWELD  software  [8]. 

Overview  of  Abaqus  Software 

Abaqus  EEA  is  a  widely  utilized  and  commercially-available,  general-purpose  EE 
software  package  developed  by  Dassault  Systemes.  The  software  enables  modeling  of  materials, 
components,  and  assemblies  under  different  loading  conditions  so  that  stress,  strain,  and  other 
relevant  quantities  can  be  determined.  Abaqus  analyses  readily  include  pre-existing 
material/structural  states  such  as  initial  residual  stresses  and  distortions.  The  reader  is  referred  to 
Abaqus  EEA  documentation  for  full  details  on  the  usage  and  input  format  for  Abaqus  software 
[9]. 


PROCEDURE 

The  transfer  of  results  from  SYSWELD  to  Abaqus  EEA  is  achieved  using  the  post¬ 
processing  tools  supplied  with  SYSWELD,  along  with  available  Abaqus  input  keywords.  A 
newly-developed  script,  SYSWELD_ABQ_field.py,  performs  this  transfer  (see  Appendix  for  a 
full  code  listing).  A  walk-through  of  the  process  is  given  below.  A  description  of  the  process  for 
generating  the  necessary  files  from  SYSWELD  is  described  followed  by  instructions  of  using  the 
script  and  including  the  translated  SYSWELD  results  in  an  Abaqus  analysis. 

The  transfer  of  SYSWELD  results  to  Abaqus  involves  several  steps  in  both  software 

tools: 

E  Generate  output  data  of  the  field  quantities  of  interest,  namely  stress  and  displacement, 
using  the  SYSWELD  post-processing  module. 

2.  Transfer  mesh  from  SYSWELD  to  Abaqus  using  SYSWELD ’s  *  .  inp  mesh  output 
capability.  Alternatively,  SYSWELD  can  read  in  the  starting  Abaqus  mesh.  It  is  critical 
that  elemental  and  nodal  numbering  be  preserved  and  that  any  mesh  generated  in  Abaqus 
be  in  the  “flat”  format;  i.e.,  not  utilize  parts  and  assemblies.  Alternatively,  SYSWELD 
analyses  can  be  performed  using  a  mesh  imported  from  Abaqus. 

3.  Translate  SYSWELD  field  quantities  to  Abaqus-interpretable  format  using 
SYSWELD_ABQ_f  ield  script. 
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4.  Run  subsequent  Abaqus  FEA  analysis  with  *  INCLUDE  keyword  to  include  initial  residual 
stress  and  distortion  input  files  generated  by  SYSWELD_ABQ_f  ield  script. 

A  full  description  of  these  steps  is  provided  below. 


SYSWELD  Results  Preparation 

The  output  required  for  the  translation  script  consists  of  a  series  of  files  tabulating  each 
stress  component  (XX,  YY,  ZZ,  XY,  XZ,  YZ)  for  every  element  in  the  form 
STRESSES_ELE_XX .  TXT,  each  distortion  component  (X,  Y,  Z)  in  the  form 
DISTORTION_NOD_X .  TXT,  and,  for  most  cases,  a  valid  exported  Abaqus  mesh.  The  steps 
described  below  document  how  to  generate  the  necessary  output  text  files  from  a  SYSWELD 
simulation  using  the  graphical  interface  of  Visual  Environment  10.5  (including  Visual  Mesh, 
Visual  Weld,  and  Visual  Viewer).  The  simulation  results  examined  were  generated  using 
SYSWELD  2015. 


SYSWELD  Mesh  Output 

The  mesh  from  the  original  CAD  file  may  have  been  modified  during  the  weld  modeling 
process  within  SYSWELD.  Therefore,  after  the  simulation,  open  the  <file  name>.vdb  file 
within  the  Visual  Mesh  program  of  SYSWELD.  To  export  the  mesh,  go  to  the  ‘Eile’  menu  and 
click  on  ‘Export’  (see  Figure  la).  This  brings  up  a  new  window  (see  Figure  lb)  enabling  the 
user  to  save  the  mesh  in  numerous  file  formats;  in  this  case,  ‘ABAQUS  files  (*.inp)’  should  be 
selected. 
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Figure  1.  Screen  shots  of  the  Visual  Mesh  function  within  SYSWELD  showing  how 
to  a)  export  mesh  and  b)  save  to  the  appropriate  Abaqus  file  format. 
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Distortion  Output 

The  distortion  of  the  mesh  after  the  SYS  WELD  simulation  can  be  viewed  by  opening  the 
<f  ile  name>_V_POST2  000  .  fdb  file  within  the  Visual  Viewer  program  of  SYSWELD. 

To  view  the  distortion,  go  to  the  ‘Results’  menu  and  click  on  ‘Contour’  (see  Figure  2).  This 
brings  up  a  new  window  (see  Figure  3).  Depending  on  settings,  the  user  must  click  differently 
to  obtain  the  same  resulting  views.  Under  the  ‘Tree’  menu,  if  set  to  Quantity,  open  the 
‘Kinematics’  drop  down  menu,  and  set  the  ‘Displacement’  drop  down  menu  to  either  X,  Y,  or  Z 
based  on  the  principle  nodal  direction  under  investigation  (see  Figure  3a).  Under  the  ‘Tree’ 
menu,  if  set  to  Entity,  open  the  ‘Node’  drop  down  menu,  and  set  the  ‘Displacement’  drop  down 
menu  to  either  X,  Y,  or  Z  based  on  the  principle  nodal  direction  under  investigation  (see  Figure 
3b). 
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Figure  2.  Screen  shot  of  the  Visual  Viewer  function  within  SYSWELD  showing  how 
to  view  the  ‘Contour’  window  used  to  display  distortion  and  stress  levels 
among  other  predicted  properties. 


The  distortion  results  can  be  output  in  text  format  by  clicking  to  the  final  step  of  the 
simulation.  Next,  go  to  the  ‘File’  menu  and  click  ‘Save  As  Ascii  Contour’  (see  Figure  4a). 
This  brings  up  a  new  window  (see  Figure  4b).  Under  the  ‘Entity’  menu,  ensure  the  text  is  set  to 
NODE.  Einally,  click  DISPLACEMENTS_NOD_X  (or  Y  or  Z),  and  save  the  output  data  to  an 
appropriate  location. 
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Figure  3.  Screen  shots  of  the  ‘Contour’  window  showing  how  to  display  the  mesh 

displaeement  depending  on  if  the  ‘Tree’  menu,  highlighted  with  dotted  red 
boxes,  is  set  to  a)  Quantity  or  b)  Entity. 
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Figure  4.  Screen  shots  of  the  Visual  Viewer  function  showing  a)  drop-down  menus  to 
save  an  ASCII  contour  (text -based  file)  and  b)  different  properties  available 
to  be  saved  as  such.  For  this  case,  DISPLACEMENTS_NOD_X  (or  Y  or  Z), 
as  highlighted  by  the  dotted  red  box,  are  required. 
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Stress  Output 

The  predicted  stress  from  the  welding  simulation  can  be  viewed  by  opening  the  <f  ile 
name>_V_POST2  0  0  0  .  f  db  file  within  the  Visual  Viewer  program  of  SYSWELD.  Similar  to 
displaying  the  nodal  distortion,  to  view  the  stress  levels  go  to  the  ‘Results’  menu  and  click  on 
‘Contour’  (as  in  Figure  2)  which  brings  up  a  new  window.  Depending  on  settings,  the  user  must 
click  differently  to  obtain  the  same  resulting  views.  Under  the  ‘Tree’  menu,  if  set  to  Quantity, 
open  the  ‘Stress’  drop  down  menu,  and  set  the  second  ‘Stress’  drop  down  menu  to  XX,  YY,  ZZ, 
XZ,  YZ,  or  XY  based  on  the  principle  stress  direction  under  investigation  (see  Figure  5a).  Under 
the  ‘Tree’  menu,  if  set  to  Entity,  open  the  ‘Solid’  drop  down  menu,  and  set  the  ‘Stress’  drop 
down  menu  to  XX,  YY,  ZZ,  XZ,  YZ,  or  XY  based  on  the  principle  stress  direction  under 
investigation  (see  Figure  5b). 


a) 


Figure  5.  Screen  shots  of  the  ‘Contour’  window  showing  how  to  display  the  predicted 
stress  depending  on  if  the  ‘Tree’  menu,  highlighted  with  dotted  red  boxes,  is 
set  to  either  a)  Quantity  or  b)  Entity. 


The  stress  results  can  be  output  in  text  format  by  clicking  to  the  final  step  of  the 
simulation.  Next,  go  to  the  ‘File’  menu  and  click  ‘Save  As  Ascii  Contour’  (see  Figure  4a); 
this  is  the  same  as  for  distortion  output.  For  stress,  however,  in  the  new  window  (see  Figure  6) 
which  opens,  ensure  the  text  is  set  to  SOLID  under  the  ‘Entity’  menu.  Finally,  click 
STRESSES_ELE_XX  (or  YY,  ZZ,  XZ,  YZ,  or  XY),  and  save  the  output  data  to  an  appropriate 
location. 
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Entity  types  Applied  on 

I  SOLID  V  On  Model  v 
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Save  Close 


Figure  6.  Screen  shot  of  the  Visual  Viewer  function  showing  the  different  properties 
available  to  be  saved  as  an  ASCII  contour  (text-based  file).  To  save  stress 
data,  STRESSES_ELE_XX  (or  YY,  ZZ,  XZ,  YZ,  or  XY),  as  highlighted  by 
the  dotted  red  box,  are  required. 


Running  Translation  Script 

Once  the  stress  and  distortion  output  is  generated  from  SYSWELD,  the 
SYSWELD_ABQ_field.  py  may  be  run.  All  SYSWELD  output  must  be  in  the  same  directory. 
The  script  itself  can  be  run  either  using  Abaqus’  built-in  Python  shell  or  any  external  installation 
of  Python  that  includes  the  csv  and  argparse  modules.  The  script  was  tested  to  work  in 
Abaqus  6.14  with  no  additional  software  beyond  a  valid  Abaqus  PEA  installation. 

The  SYSWELD_ABQ_field.py  script  has  separate  functionalities  for  nodal  and  element 
integration  point  quantities  and  must  be  executed  twice:  once  for  stress  and  once  for  distortion. 
The  stress  and  distortion  quantities  are  handled  as  follows: 

•  Residual  stresses,  described  at  element  integration  points,  are  processed  and  formatted 
using  the  *INITIAL  CONDITIONS ,  TYPE=STRESS  Abaqus  keyword.  SYSWELD 
generates  an  individual  file  for  each  component  of  stress  (<Txx’  <^yy>  <^zz>  ^xy>  ^xz>  ^yz)- 
The  script  generates  a  standalone  Abaqus  input  file,  stress .  inp,  that  provides  the 
individual  stress  components  in  a  six-column  vector  that  is  tabulated  for  each  element. 
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•  Initial  distortions  or  displacements,  a  nodal  quantity  with  three  components  (X,  Y,  Z)  are 
converted  to  the  *  imperfection,  input=filename  Abaqus  keyword  format. 
SYSWELD  generated  an  individual  file  for  each  component  of  displacement 
(Ax,  Ay,  Az) .  A  file  displacement .  csv  is  generated  that  assembles  these  individual 
displacement  components  into  a  three-column  vector  associated  with  each  node  number. 

The  SYSWELD_ABQ_field.py  script  is  invoked  through  a  command-line  driven  interface  using 
the  syntax: 

abaqus  python  SYSWELD_ABQ_field. py  -fieldType  <displacement | stress> 
-path  <path  to  SYSWELD  field  data>  -scaleFactor  <SF> 

The  arguments  describe: 

•  fieldType:  Field  of  interest.  Stresses  are  six-value  tensor  quantities  at  each  element 
whereas  initial  displacements  are  three-value  vectors  at  each  node. 

•  path:  Path  to  Sysweld-generated  stress  or  displacement  data. 

•  scaleFactor:  Scale  factor  to  convert  SYSWELD  units  to  units  in  Abaqus  analysis 
for  field  of  interest. 

A  help  description  can  be  obtained  by  executing: 

abaqus  python  SYSWELD_ABQ_field. py  -h 

The  commands: 

python  SYSWELD_ABQ_field.py  -fieldType  stress  -path 
/myDirectory/results  -scaleFactor  le6 

python  SYSWELD_ABQ_field.py  -fieldType  displacement  -path 
/myDirectory/results  -scaleFactor  0.001 

will  generate  stress .  inp  and  displacement .  csv  files  using  the  SYSWELD  stress  and 
displacement  data  located  in  /myDirectory/results.  A  multiplier  of  le6  and  0.001  will  be 
applied  to  all  stress  values  and  displacement  values,  respectively.  It  is  the  user’s  responsibility  to 
ensure  that  consistent  unit  systems  are  maintained.  If  a  component  is  not  defined  by  the  data 
provided  from  SYSWELD,  a  null  value  is  assumed  to  ensure  either  a  three-column  vector  for 
displacement  or  a  six-column  stress  vector  for  stress  is  obtained. 


It  is  strongly  recommended  that  a  static  step  be  undertaken  with  the  SYSWELD- 
generated  residual  stresses  and  distortions  to  ensure  precise  equilibrium.  Once  this  step  is  taken, 
the  Abaqus  analysis  of  interest  can  be  performed. 

USAGE  EXAMPLE 

The  transfer  of  SYSWELD  results  to  Abaqus  was  demonstrated  for  a  B2V.1  joint,  which 
is  a  single-V  butt  joint  without  a  backing  bar  [10].  The  simulated  material  was  A1  5052  using  the 
gas-metal  arc  welding  (GMAW)  process.  A  simulation  of  the  welding  process  of  this  joint  was 
conducted  in  SYSWELD  and  is  documented  in  another  NSWCCD  report  [11].  While  this 
specific  weldment  does  not  represent  an  actual  shipboard  application,  it  was  selected  because  it 

13 


NSWCCD-61-TR-20 15/31 


demonstrated  numerous  weld  processing  conditions,  including  multiple  weld  passes,  flipping  the 
weldment,  and  a  back-gouge,  in  order  to  validate  the  SYSWELD  to  Abaqus  technique.  Figure  7 
shows  a  view  of  the  mesh,  clamping  conditions,  and  the  three  weld  passes,  shown  in  purple, 
pink,  and  gold  coloring. 


Figure  7.  SYSWELD-generated  model  of  an  aluminum  B2V.1  joint  showing  the 
initial  mesh  and  clamping  conditions. 


Translating  Results 

Once  the  SYSWEED  output  steps  are  complete,  the  SYSWELD_ABQ_field.py  script  is 
run  to  generate  a  stress. inp  file.  This  file  is  directly  included  in  the  Abaqus  analysis  including  the 
following  keyword  in  the  main  input  file: 

*Include ,  input=stress . inp 

A  contour  plot  showing  the  final  residual  stress  state  from  SYSWEED  simulations  and 
the  initial  state  of  an  Abaqus  analysis  using  these  stresses  is  show  in  Figure  8.  Note  that  the 
mesh  was  generated  in  SYSWEED,  exported  to  Abaqus,  and  scaled  from  mm-kg-MPa  to  MKS 
units.  The  stresses  were  translated  accordingly.  Prior  to  undertaking  further  Abaqus  analysis,  it  is 
strongly  recommended  that  a  static  equilibrium  step  be  performed  to  ensure  exact  static 
equilibrium. 
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Abaqus  FEA  Input 


Figure  8.  Comparison  of  Von  Mises  stress  levels  from  a)  SYSWELD  simulation 
output  values  and  b)  same  values  translated  to  Abaqus  FEA.  Minor  color 
variations  are  due  to  differences  in  contour  coloring. 

Using  Initialized  Local  Weld  Model  in  Structural  FEA 

With  the  ability  to  include  initial  residual  stress  and  distortion  from  SYSWEED 
simulations  in  Abaqus  structural  analyses,  global  loads  can  readily  be  transferred  into  the  local 
weld  detail  using  the  sub-modeling  capability  available  in  Abaqus  FEA.  This  capability 
automatically  maps  the  edge  displacements  of  the  coarse  global  mode  to  edge  boundary 
conditions  on  the  local  model. 

An  example  application  of  this  is  shown  in  Figure  9.  The  initial  state  of  the  local  weld 
and  the  unloaded  state  are  shown  in  Figure  9a.  As  displacement  is  applied  to  the  edges  of  the 
global  structural  model,  a  complex  stress  state  in  the  weld  is  observed  to  develop  in  Figure  9b. 
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This  result  is  very  different  from  the  result  obtained  without  including  residual  stresses  from 
welding  due  to  the  high  level  of  these  stresses  relative  to  the  stress  induced  by  structural  loading. 


Initial  Residual  Stresses 


Loaded  Stress  State 


Figure  9.  Embedded  weld  with  initial  residual  stresses,  determined  by  SYSWELD,  in 
global  shell  mesh  a)  prior  to  a  loading  simulation  in  Abaqus,  and  b)  post¬ 
loading  in  Abaqus.  Notice  the  stresses  in  the  weld  and  global  shell  mesh 
under  uniaxial  tension  loading  along  stiffener  direction.  Coupling  from 
global  shell  mesh  to  local  weld  mesh  is  fully  automated  using  one-way. 

CONCLUSIONS 

A  process  and  necessary  software  for  translating  SYSWEED  results  from  inclusion  in 
Abaqus  EEA  structural  analyses  has  been  developed  and  its  usage  demonstrated  the  successful 
integration  of  the  two  software  programs.  The  application  of  this  newly-developed  capability  to 
structural  analysis  of  a  multi-pass  aluminum  weld  illustrates  the  inability  of  structural  EEA 
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calculations  to  capture  local  weld  detail  response.  As  failures  often  initiate  at  these  locations,  the 
U.S.  Navy  is  now  able  to  investigate  how  the  welding  process  modifies  the  resultant  fracture  and 
fatigue  behavior  for  structural  components  on  its  vessels. 
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APPENDIX:  SYSWELD_ABQ_field.py  CODE 

#  SYSWELD_ABQ_field.py 

#  Written  by  Ken  Nahshon  {k.en.nahshon@navy.mil).  Naval  Surface  Warfare  Center 

#  Carderock  Division 

#  27  July  2015 

#  This  code  is  subject  to  all  USC  regulations.  NSWCCD  shall  not  be  held  liable  for  any  use  or 

#  misuse  of  this  code.  The  code  is  stictly  provided  as-is. 

#  - 

#  Converts  field  quantities  from  SYSWELD  to  ABAQUS 

#  Input  consists  of  a  directory  with  SYSWELD  results  e.g.  DISPLACEMENTS_NOD_X.txt  etc  and  STRESSES_ELE_XX.txt 

#  Output  consists  of  a  file  stress. inp  consisting  of: 

# 

#  stress. inp:  *INITIAL  CONDITION  pre-stress  values.  The  stress  file 

#  can  be  included  using  *Include,  input=stress . inp 

#  or  directly  copied  and  pasted  into  the  ABAQUS  input  file. 

# 

#  displacement . CSV :  delta  x,  y,  and  z  for  moving  nodal  coordinates. 

#  ABAQUS  can  move  the  nodal  coordinates 

#  using  *IMPEREECTION  input=displacement . CSV 

#  - 

#  Usage: 

# 

#  python  SYSWELD_ABQ_f ield.py  -fieldType  <field_type>  -path  <path  to  SYSWELD  files>  -scaleEactor  <scaleEactor> 

#  Where, 

#  field_type  can  be  either  'stress'  or  'displacement'.  Currently,  other  fields  are  not  supported 

#  scaleEactor  is  a  floating  point  value  multiplying  the  SYSWELD  results  to  account  for  unit  changes 

#  working_directory  is  the  location  of  SYSWELD  output  files 

# 

#  Example: 

# 

#  python  SYSWELD_ABQ_f ield.py  -fieldType  stress  -scaleEactor  le6  -path  /runs/weldcalc/ 

# 

#  Will  execute  the  script  to  generate  a  stress. inp  file  in  /runs/weldcalc/  using  the  SYSWELD  output  located 

#  in  that  directory  with  stresses  from  SYSWELD  multiplied  by  the  scale  factor  le6 

# 

# - 

import  os 
import  CSV 
import  sys 
import  argparse 

def  f ilesWithString { f ileStr,  directory): 

"""Returns  a  list  of  filenames  with  fileStr  in  name  in  directory\n 
Usage:  f ilesWithString (fileStr,  directory) """ 
listing=os . listdir (directory) 

return  [elem  for  elem  in  listing  if  elem. find { fileStr) >0 ] 

def  openInputEile { inputEileName,  readWriteOption) : 

"""Safe  file  open  with  exceptions""" 
import  sys 

print  'Processing:  ',  inputEileName 
try: 

inputEile=open (inputEileName,  readWriteOption) 
return  inputEile 
except  lOError  as  e: 

print  "Error  in  file  open  in  openInputEile : \n" 
print  "I/O  error{{0}):  { 1 }". format (e . errno,  e.strerror) 
except : 

print  "Unexpected  error:",  sys . exc_inf o { ) [ 0] 
raise 

def  openCsv { f ileName ,  commentChar  =  {'#',  '%')/  **kwargs) : 

Read  a  .csv  file  with  fileName  (string)  and  return  a  list  of  lines  where  each  line  is  a  list 
of  entries  seperated  by  commas  in  fileName.  Routine  strips  out  lines  beginning  with 
commentChars  and  additional  whitespace.  Additional  comment  markers  can  be  included. 

Usage:  \n\n 

parseCSVtoList (fileName,  commentChar  =  {'#',  '!'/  '%')) 

whitespace= { '  ',  '\t')  #  Define  whitespace  to  ignore 

f=openInputEile { fileName ,  'r') 

reader=csv. reader { f ,  **kwargs) 

lineCount=0 

inputData=  [  ] 

for  row  in  reader: 

if  {len{row)  ==  0)  or  {row[0] [0]  in  commentChar) :  #  Check  for  comment  characters  and  skip  entire  line 
continue  #  Skip  empty  lines 

else : 
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#  Eliminate  additional  whitespaces 
for  char  in  whitespace: 

row= [row [index] . replace (char,  "")  for  index  in  range{0,  len{row))] 
inputData. append (row) 
lineCount  +=  1 
return  inputData 


#  MAIN  CODE 

#  ********** 


#  Show  title  info 


titleBlock  =  ""'i\j^********************************************************************** 
\nSYSWELD_ABQ_f ield.py,  a  tool  to  convert  SYSWELD  fields  to  ABAQUS  formats. 
Currently  supports  stress  and  displacment  types. 


#print  titleBlock 

parser  =  argparse.ArgumentParser {description=titleBlock) 

parser . add_argument { ' -f ieldType ' ,  required=True,  metavar= { ' < displacement  I stress> ' ) ,  type=str ,  nargs=l, 
help= ' Specif y  field  type') 

parser . add_argument {' -path ' ,  required=True,  metavar= { ' <Path  to  SYSWELD  field  data>'),  type=str,  nargs=l, 
help= ' Specif y  directory  for  SYSWELD  files') 

parser . add_argument { ' -seal eEac tor ' ,  required=True,  metavar= { ' <SE> ' ) ,  type=f loat,  nargs=l ,  help= ' Specify  scale 
factor  for  field') 

parser ._optionals . title  =  "Elag  arguments" 

#  Parse  arguments 

args  =  parser. parse_args {) 

#  Get  command  line  input 
fieldType  =  args . f ieldType [0 ] 
scaleEactor  =  args . scaleEactor [0 ] 
directory  =  args.path[0] 

#  Make  sure  directory  exists 

if  os. path. exists (directory)  is  Ealse: 

print  Eatal  Error,  %s  is  not  a  valid  path  ***'  %directory 

sys . exit (0 ) 

#  Make  sure  /  is  last  character  in  directory  name 

#if  directory [-1]  '/': 

#  directory  =  directory  +  ' / ' 

#  Precalculation  output 

print  '\nEield  type:  %s'  %fieldType 

print  'Multiplying  SYSWELD  %s  by  scale  factor:  %s'  % (fieldType,  scaleEactor) 
print  'Working  directory:  %s\n'  %directory 

#  Switch  to  working  directory 
os . chdir (directory) 

#  - 

#  Set  up  valid  field  types,  directions 

#  - 

# 

#  Valid  field  types 

fieldList  =  ['displacement',  'stress'] 

#  List  of  directions  for  tensors  (elements) 

directionList= [ 'XX ' ,  ' YY ' ,  'ZZ',  'XY',  ' YZ ' ,  ' ZX ' ] 

#  List  of  directions  for  nodes 

nodalDirectionList= [ 'X ' ,  'Y',  'Z'] 

#  Check  for  supported  field  type,  set  file  prefix  for  SYSWELD.  If  not  supported,  exit, 
if  fieldType  ==  'stress': 

prefix  =  ' STRESSES_ELE_' 
elif  fieldType  ==  'displacement': 

prefix  =  'DISPLACEMENTS_NOD_' 
else : 

print  '**  ERROR:  following  field  type  is  not  supported:  ',  fieldType 
sys . exit ( ) 

#  Perform  actions  based  on  field  type 

#  Stress  translation 

if  fieldType  ==  'stress': 

#  Set  a  dictionary  of  stress  values,  read  stresses  from  direction  if  file  exists 
stress  =  dictO 

for  direction  in  directionList : 

fileName  =  prefix  +  direction  +  '.txt' 
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if  fileName  in  os . listdir { directory) : 

data=openCsv (directory+f ileName,  delimiter= ' \t ' ,  skipinitialspace  =  'True') 

#  Strip  header  off  data  file.  Currently,  SYSWELD  fills  first  8  lines 
data^data [ 8 : ] 

#  Read  element  list 

elList  =  [  int (dataLine [0] )  for  dataLine  in  data] 

#  Read  stresses,  break  into  directions  using  dictionary 

stress [direction]  =  [ float (dataLine [1 ]) *scaleFactor  for  dataLine  in  data] 

numElements  =  len{elList) 

#  Backfill  with  zeroes  for  stress  components  that  don't  exist 
for  direction  in  directionList : 

if  direction  not  in  stress . keys {) : 

stress [direction]  =  [0.0  for  line  in  range{0,  numElements)] 

#  write  out  data  as  six  column  vector 

with  open (directory  +  ' stress . inp ' ,  'wt')  as  f: 

f . write ('* *Autogenerated  prestress  from  SYSWELD\n') 
f . write ('* Initial  Conditions,  type='  +  fieldType  +  '\n') 
for  i  in  range ( 0 , numElements ) : 

line-  '{0},  {1},  {2},  {3},  {4},  {5},  { 6 } \n '. format ( str ( elList [ i ]) , 

St r ( stress [' XX ' ] [i] ) ,  str ( stress [' YY ' ] [i] ) ,  str ( stress [' ZZ ' ] [i] ) , 

str ( stress [' XY ' ]  [i] ) ,  str ( stress  [' YZ ' ]  [i ] ) ,  str ( stress [' ZX ' ]  [i] ) ) 
f . write ( line ) 

print  '\nWrote  %s  in  directory  %s\n'  %(' stress . inp ' ,  directory) 

if  fieldType  ==  'displacement': 
displacement  =  dictO 
nodeList  =  listO 

for  direction  in  nodalDirectionList : 

fileName  =  prefix  +  direction  +  '.txt' 
if  fileName  in  os . listdir ( directory) : 

data=openCsv (directory+f ileName,  delimiter- ' \t ' ,  skipinitialspace  =  'True') 

#  Strip  header  off  data  file 
data-data [ 8 : ] 

nodeList  =  [  int (dataLine [0] )  for  dataLine  in  data] 

displacement [direction]  =  [ float (dataLine [1 ]) *scaleEactor  for  dataLine  in  data] 

numNodes  =  len (nodeList) 

#  Backfill  with  zeroes 

for  direction  in  nodalDirectionList: 

if  direction  not  in  displacement . keys () : 

displacement [direction]  =  [0.0  for  line  in  range(0,  numNodes)] 

#  write  out  data 

with  open ( directory  +  ' displacement . csv ' ,  'wt')  as  f: 

f . write ('* *Autogenerated  initial  distortions  from  SYSWELD\n') 

#  f . write (' *Initial  Conditions,  type-'  +  fieldType  +  '\n') 

for  i  in  range ( 0 , numNodes) : 

line  -  '{0},  {1},  {2},  { 3 } \n '. format ( str (nodeList [ i] ) ,  str (displacement [ 'X '] [i] ) , 
str (displacement [ ' Y ' ] [ i] ) ,  str (displacement [ ' Z ' ] [i] ) ) 
f . write ( line ) 

print  '\nWrote  %s  in  directory  %s\n'  %(' displacement . csv ' ,  directory) 
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